1. Joint Modeling
1.1 Trace plots for convergence check
The current MCMC setting is:
- 100,000 iteration;
- 90,000 burn-in;
- 10 thinning.
1.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.00 1.00
## LDevsum 1.01 1.03
## dh0 1.25 1.68
## dh1 1.25 1.68
## dh2 1.00 1.01
## dl0 1.06 1.18
## dl1 1.18 1.47
##
## Multivariate psrf
##
## 1.18
1.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
1.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1141.77712 | 22552.0905 |
| DIC3 | 1111.55453 | 22638.2755 |
| PWAIC | 22.66121 | 203.6814 |
| WAIC | 1122.82311 | 22651.9122 |
2. Separate Modeling
2.1 Trace plots for convergence check
The current MCMC setting is:
- 100,000 iteration;
- 90,000 burn-in;
- 10 thinning.
2.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.01 1.04
## LDevsum 1.00 1.01
## dh0 5.56 11.07
## dh1 10.17 18.41
## dh2 5.84 10.29
## dl0 1.11 1.32
## dl1 1.02 1.05
##
## Multivariate psrf
##
## 9.02
2.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
2.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1388.1021 | 33774.6195 |
| DIC3 | 1244.4182 | 33891.6284 |
| PWAIC | 89.7843 | 169.8478 |
| WAIC | 1312.5302 | 33896.8859 |